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STABILITY  OF  TWO-DIMENSIONAL  RECURSIVE  FILTERS 


0 


$ 


I. 


Introduction 


In  digital  image  processing,  we  often  want  to  do  linear  filtering.  The  use  of 

recursive  instead  of  non-recursive  filters  has  the  potential  of  saving  computation  time. 

A  two-dimensional  digital  recursive  filter  is  characterized  by  the  two-dimensional 

z-transform  R 

2  a  z  z0 

™  -n  «=n  mn  1  ^ 

-  (1) 
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where  a  and  b  are  constants,  and  without  loss  of  generality  we  can  set  b  =  1. 
mn  mn  °  y  oo 

The  degrees  of  the  numerator  and  the  denominator  polynomials  do  not  have  to  be  equal, 

since  some  of  the  coefficients  a  and  b  can  be  zero.  The  variables  z,  and  z~ 

mn  mn  1  l 

are  defined  as 
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-s2B 


(2) 


where  s^  and  s,;  are  respectively  the  horizontal  and  vertical  complex  spatial  frequency 
variables  and  A  and  B  are  constants  (sampling  periods  in  the  horizontal  and  vertical 
directions,  respectively). 

Let  f(m,n)  andg(m,n)  be  the  input  and  output,  respectively  of  the  filter.  Then 
the  spatial -domain  difference  equation  corresponding  to  Eq.  (1)  is 
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Fig.  1.  Four  directions  of  recursion. 
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Fig.  2.  Initial  conditions  for  a  casual 
filter. 
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2  2  b  g  (M-m,  N-n)  =  2 

m=0  n=0  m=0 

We  can  express  g  (M,N)  in  terms  of  the  rest,  and  thus  obtain  a  recursive  filter  which 

recurses  in  the  (Em,En)  direction.  We  can  also  express  g  (M-p,N-q)  in  terms  of  the 

rest,  and  obtain  a  recursive  filter  which  recurses  in  the  (-m,  -n)  direction,  Similarly, 

we  can  get  from  Eq.  (3)  recursive  filters  recursing  in  the  (Em,  -n)  and  the  (-m,+n) 

directions.  See  Fig.  1.  We  shall  call  the  recursive  filter  that  recurses  in  the  (Em, 

En)  direction  causal.  Thus,  the  z-transform  of  Eq.  (1)  can  be  associated  with  four 

different  recursive  filters.  However,  unless  otherwise  specified,  we  shall  associate 

it  with  the  causal  one. 

To  start  recursing,  we  need  initial  conditions.  For  the  causal  filter,  we  need 
to  know  the  values  of  the  output  g  in  the  shaded  region  in  Fig.  2. 

There  are  two  major  problems  in  the  design  of  a  recursive  filter:  approximation 
and  stability.  The  approximation  problem  consists  of  determining  the  coefficients 
amn  anci  ^mn  in  Eq.  (1)  so  that  H(e  e  approximates  a  given  frequency  response 

H^(ju,  jv),  where  u  and  v  are  respectively  the  horizontal  and  vertical  spatial  frequencies. 
In  order  for  the  resulting  filter  H(z^,  z^)  to  be  useful,  we  have  to  require  that  it  be 
stable,  i.e.  ,  if  we  expand  H(z^,  z^)  into  a  power  series 

CO  00 

H(z  ,  z„)  =  2  2  h  zm  z9  (4) 

1  4  n  mn  12 


q 

2  a  f(M-m,N-n)  (3) 
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the  coefficients,  h^,  which  is  the  impulse  response  of  the  filter,  should  satisfy 


oo  a 

2  2 

m=0  n=0 


h  <oo 

mn  1 


(5) 


If  the  filter  is  unstable,  any  noise  (including  the  roundoff  errors  in  computation)  will 
propagate  through  the  output  and  possibly  be  amplified.  Furthermore,  in  many  important 
applications,  such  as  inverse  filtering  for  image  restoration,  the  initial  conditions 
are  usually  not  known.  We,  therefore,  desire  that  the  result  of  the  recursive  filtering 
be  roughly  independent  of  the  initial  conditions.  This  will  be  true  if  the  impulse 
response  I  h  I  wO  for  m,  n  >  A,  where  A  is  a  constant  which  is  much  smaller 
than  the  linear  dimensions  of  the  input  and  output  images.  Therefore,  in  some 
cases,  stability  itself  is  not  sufficient;  we  require  in  addition  that  the  impulse  response 
die  out  fast  enough. 

In  this  paper,  we  shall  discuss  some  aspects  of  the  stability  problem  of  two- 

dimensional  recursive  filters.  In  particular,  we  shall  derive  a  simplified  version 

1  2 

of  Shanks  stability  theorem  ’  ,  and  show  that  it  is  equivalent  to  Ansell  stability 

..  3,4 

theorem 

II.  Shanks  Stability  Theorem 

Both  the  theorems  of  Shanks  and  Ansell  tell  us  stability  conditions  in  the  fre¬ 
quency  domain.  In  this  section,  we  shall  state  and  prove  Shanks  stability  theorem 
and  discuss  its  implications. 
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a)  Statement  of  the  Theorem 


Theorem  1  (Shanks) 

A  causal  recursive  filter  with  the  z-transform  H(Zj,  z ,^)  =  A(z^,  z^/Bfc^.z.,), 
where  A  and  B  are  polynomials  in  z^  and  Z£.  is  stable,  if  and  only  if  there  are  no 
values  of  z^  and  Z2  such  that  B(z^,  Z2)  =  0  and  |  z^  j  £  land  |  z ^  |  ^  1. 
b)  Proof 


Let 


H(zlt  z2)  = 


A(z1(  z2) 

B(z1>  z2) 


CO  CO 

=  2  2 
m=0  n=0 
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mn 
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z2 


(6) 


00  CO 

What  we  want  to  show  is  that  2  2  |  h  |  <  00  if  and  only  if  H(z.,  z9)  is  analytic 

m=0  n=0  mn 

in  the  region  D  =  4 (z ^ ,  z2);  |  Zj  |  £  1  n  |  z2  |  ^  1 

The  "if"  part:  If  H(z^,  Z2)  is  analytic  in  D,  we  can  find  « >  0  such  that 
H(z  j,  x^)  is  analytic  in  D  ^  =  j(z  ^ ,  x^)  '•  |z^  |<l+en|z2|<'+e^,  which 

implies  that 

2  2  h  x.  Xn,  is  absolutely  convergent  in  D,. 

m  n  mn  12  3  b  1 

Therefore  2  2  I  h  I  <». 

m  n  1  mn  1 

The  "only  if"  part:  If  2  2  I  h  I  <  00,  then  by  the  M-test, 

3  r  m  n  1  mn  1  3 


2^  2  h  x™  x 2  is  absolutely  convergent  in  D,  which  implies  in  turn  that  H(z^,Z2) 
is  analytic  in  D.  Q.  E.D. 


¥ 
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c)  Test  Procedure 


To  test  stability  using  Shanks  theorem  is  hard  work.  One  way,  e.g.  ,  is  to 
map  the  unit  disk  d^  =(  z^  |  Zj  |  £  1)  in  the  z^  -plane  into  the  z2  -plane  by  the  implicit 
mapping  relation  B(z^,  z2)  =  0.  The  filter  is  stable  if  and  only  if  the  image  of  d^  in 
the  z2-plane  does  not  overlap  the  unit  disk  d2  =  I  s  *)•  This  is  hard  work, 
because  for  each  particular  value  of  z^  =  z°  ed^,  we  have  to  solve  the  equation 
B(Zj  ,  z2)  =  0  for  z2-  And  we  have  to  do  that  for  all  (in  practice,  a  large  number 
of  points) in  d 

d)  An  Implication 

We  emphasize  that  Shanks  theorem  was  stated  for  causal  filters,  that  is, 
filters  which  recurse  in  the  (+m,+n)  direction.  Let  us  consider  the  filter  associated 
with  H(Zj,  z2)  =  A(Zj,  z2)/B(Zj,z2)  which  recurses  in  the  (-m, -n)direction.  We 
assume  that  the  orders  of  A  in  both  z^  and  z2  are  smaller  than  the  corresponding 
orders  of  B.  Then,  the  stability  condition  becomes:  The  image  of(Zj,|  zj^  1)  in 
the  z2-plane,  according  to  B(z^,  z2)  =  0,  does  not  overlap  (z2.  |z2  |  s  1.) 

This  can  be  shown  by  expanding  H(z^,  z2)  in  terms  of  z^  and  z^\ 


A(zj,  z2) 

H(zr  z2)  z2) 


oo 

=  2 


00 
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m=0  n=0 


mn 


-m  -n 
zl  z2 
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Making  the  change  of  variables 

-1 

W1  =  Z1 

-1 


w2  =  z2 


(8) 
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\- 
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We  get 


m  n 

i  w0 


(9) 


m  n  mn 


Now 


A1  (wi>  w2^ 


B1  (w  x ,  w2) 


(10) 


Where  a  and  b  are  non-negative  integers,  and  =  0  if  and  only  if 

B(w^,  W2  S  =  0.  Applying  Shanks  theorem  to 


(ii) 


Bi(wr  w2) 


we  know  that  2  2  I  h*  I  <  oo  if  and  only  if  the  image  of  (w.,  I  w,  I  £  1)  in  the 

m  n  1  mn  1  J  b  '  1  1  1  1  ' 


w 2  -plane,  according  to  B^(w^,  w^)  =  0,  does  not  overlap  (w7;  |  ^ ^  |  £  1).  But 
W1  ~Z1^’  anc^  w2  =  z2  *  our  asser^on  has  been  proved. 

Similarly,  we  can  get  the  stability  conditions  for  the  filters  associated  with 
H(z^,  z^)  that  recurse  in  the  (+m,  -n)  and  the  (-m,-Hn)  directions.  We  can  easily 
convince  ourselves  that  at  most  one  of  the  four  stability  conditions,  corresponding 
to  the  four  recursing  directions,  can  be  satisfied*  Therefore,  we  have  the  following 
Theorem  2 

Among  the  four  recursive  filters  we  can  associate  with  the  z-transform  of 
Eq.  (1),  or  equivalently  the  difference  equation  of  Eq.  (3),  at  most  one  can  be  stable. 
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SINGULAR  z,-  PLANE  B^.z^O 

POINTS  OF  or 


z2- PLANE 


Fig.  3.  The  mapping  z ^  -  f(Zj). 
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III.  A  Simplified  Version  of  Shanks  Theorem 

Shanks'  result  as  stated  in  Theorem  1  is  very  tedious  to  apply.  We  now 
show  that  his  stability  conditions  can  be  simplified  considerably. 

a)  The  Result 
Theorem  3 

A  causal  filter  with  a  z-transform  H(z^,  z^)  =  A(Zj,  z0)/B(Zj,  z^), 
where  A  and  B  are  polynomials,  is  stable,  if  and  only  if:  (1)  the  map  of  £d^  e 
(Zy  |zjJ  =  1)  in  the  Z2  “plane,  according  to  B(z^,  z9)  =  0, lies  outside  d2  z  (z2*>  \7 2 

£  1);  and  (2)  no  point  in  d^  =  (z^;  |z^  |  £  1)  maps  into  the  point  z9  =  0  by  the  relation 
B(Zj,  z2)  =  0. 

b)  Proof 

We  want  to  establish  that  the  stability  conditions  of  Theorems  1  and  3  are 
equivalent.  That  the  stability  conditions  of  Theorem  1  imply  those  of  Theorem  3 
is  obvious.  So  we  proceed  to  show  the  implication  in  the  other  direction. 

The  two-variable  polynomial  B(z^,  Z2)  =  0  defines  an  algebraic  function‘d 
z2  =  ^zi^*  We  ^rst  modify  the  unit-circle  contour  in  the  z^  -plane  to  exclude  any 
singular  points  of  f  inside  the  contour,  resulting  in  a  modified  contour  0dj  as  shown 
in  Fig.  3.  We  use  d^  to  denote  the  closed  region  enclosed  by  0d^.  A  point  z^  =  z  ^ 
is  called  a  singular  point  of  Z2  =  f(Zj),  if  B(z^,  Z2)  =  0,  considered  as  an  equation 
in  z9,  has  multiple  (finite  or  infinite)  roots. 
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According  to  the  theory  of  algebraic  functions,  in  dj,  the  function  z 2  =  f(z^) 
has  a  number  of  branches,  each  of  which  is  holomorphic.  Therefore,  from  the  max¬ 
imum-modulus  theorem,  the  maximum  of  |f(z^)|  over  d^  occurs  on8d^,  and  the  min¬ 
imum  of  |f(Zj)|  over  dj  can  occur  in  the  interior  only  if  the  minimum  is  zero.  How¬ 
ever,  condition  (2)  of  Theorem  3  says  f(z^)  is  never  zero  in  d^.  Therefore,  the  min¬ 
imum  of  f(Zj)  occurs  on  8dj,  i.e.  , 

If  (dp  |  :>  min.  (f(8dp|  (12) 

which  implies  that:  if  |f(8dp|  >1,  then  | f (dp |  >  1;  i.e.,  to  ensure  that  f (dp  lies 
outside  the  unit  disk  d2  =  (z^,  |  z2|  £  1),  it  is  sufficient  to  ensure  that  f(8d’1> 
lies  outside  d2< 

We  are  almost  there,  but  not  quite.  What  we  really  want  to  show  is  that: 
if  |f(8d  j)  |  >  1,  then  | f (d  ^ )  |  >  1,  where  d  ^  =  (z  |z  ^  |  £  1)  and  8d^  5  (z.‘,  |  z  jJ  =1). 
Since  the  detour  in  8dj  can  be  any  path  leading  from  8d^  to  the  singular  point,  what 
is  left  to  show  is  simply  that  |f(s)|  >  1  where  s  is  the  singular  point.  But  since  each 
branch  of  z2  =  f(zp  is  continuous  at  z^  =  s,  and  since  |f(s+e  e  |  >1  for  arbitrary 
small  e  and  any  0,  we  have  |f(s)  |  >  1.  Q.  E.  D. 

c)  Test  Procedure 

To  test  stability  using  Theorem  3,  we  have  to  map  8d^  =  (z  ',  |z^|  =  1)  into  the 
z2  -plane  according  to  B(z^,  z2)  =  0,  and  see  whether  the  image  lies  outside  d2  e. 

(z2i  |z2|  £  1).  Also,  we  have  to  solve  B(z^,0)  =  0  to  see  whether  there  is  any  root 
with  magnitude  smaller  than  1. 
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IV. 


Ansell  Stability  Theorem 

To  test  the  stability  of  a  two-dimensional  recursive  filter  using  Theorem  3 


is  much  simpler  than  using  Shanks’  original  Theorem  1.  However,  the  procedure 
is  still  infinite  in  the  sense  that  in  principle  we  have  to  map  every  point  on  the  con¬ 
tinuum  <9d,  into  the  z()  -plane.  In  this  section,  we  show  that  Theorem  3  can  be 
reduced  to  a  result  due  to  Ansell.  We  shall  see  that  Ansell’s  result  enables  us  to 
test  stability  in  a  finite  number  of  steps  -  which,  unforntunately,  can  still  be  very 
tedious. 


and 


By  making  the  change  of  variables: 
1  -  z, 


Pi  = 


1 


1  1 


T 


1  -  z, 


p2  i 


+  zr 


and  let 


E  (p..  P2) 

H(z  ,  Z„)  =-rr; - 

1  2  F  (pr  p2) 


(13a) 


(13b) 


(14) 


where  E  and  F  are  polynomials  in  p^  and  pr  We  can  restate  Theorem  3  as: 

Theorem  4  (Ansell) 

The  causal  recursive  filter  H(z^,  z^)  is  stable,  if  and  only  if:  (1)  for 
all  real  finite  w,  the  complex  polynomial  in  p  ,,  F(jw,  p}),  has  no  zeros  in  Re  p}  ^0; 
and  (2)  the  real  polynomial  in  p^  F(p^,  1),  has  no  zeros  in  Re  p^  £  0. 

3 

This  theorem  is  essentially  the  same  as  Lemma  6  of  Ansell  .  The  nice  thing 
about  this  theorem  is  that  condition  (1)  can  be  tested  using  standard  techniques  of 
circuit  theory*5. 
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Theorem  5  (Ansell) 


Condition  (1)  of  Theorem  4  is  equivalent  to  the  following:  Let  us  express 
F(jw,  jft),  where  w  and  ft  are  real,  as 


F(jw,  jft)  =  b  (w)  ftn  +  b.  (w)  ftn  *  +  ...+  b  (w) 

o  i  n 


+  j  [aQ(w)  ftn  +  a.l  (w)  ftn  1  +  .  . .  +  an(w)] 


(15) 


where  a.(w)  and  b.(w)  are  real  polynomials  in  w  and  either  aQ(w)  or  bQ(w)  is  not 

identically  zero.  Let  C  (w)  be  defined  as 

r ,  s 


C  =  a  b  -  a  b 
r,  s  r  s  s  r 


(16) 


for  Osr,ssn  (using  zero  for  a's  and  b's  not  present  inF(jw,jft)).  And  let  D(w)  denote 

the  n  x  n  symmetrical  polynomial  matrix  whose  typical  element  D^(w)  (1  £  i,  j  sn) 

is  the  sum  of  all  those  C  (w)  (0  sr,s  £  n)  for  which  both 

r  ,  s 


s  +  r  =  i  +  j  -  1 


(17a) 


and 

s  -  r  >  |  i  -  j  |  (17b) 

are  satisfied.  Then,  the  n  successive  principal  minors  of  D(w)  must  be  positive  for 
all  real  w. 

3  7  6 

This  is  part  of  Theorem  II  of  Ansell  ’  .  We  note  that  Sturm's  method  can 
be  used  to  test  whether  each  minor  of  D(w)  is  positive  for  all  real  w. 
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V. 


Examples 

We  now  give  some  examples  illustrating  the  use  of  Theorem  3  and  Theorems 


w' 


4  and  5.  Ideally,  we  would  like  to  determine  what  relations  the  coefficients  a.(  and 
b.j  in  Eq.  (1)  must  satisfy  to  ensure  the  stability  of  a  causal  recursive  filter  H(z^,Z2). 
Unfortunately,  this  does  not  seem  possible  except  for  very  simple  filters, 
a)  First-order  Filters 
Let  us  consider  the  first-order  filter 


H(Zj.,  z2^  -  1  +  az^  +  bz2 

where  a  and  b  are  constants.  In  this  simple  case,  we  can  easily  determine  the 
stability  conditions  in  terms  of  a  and  b  using  Theorem  1. 

Let 

B(z^,  Z2)  =  1  +  az^  +  bz2  =  0 
we  have 

1  a 

z2  “  "  b  "  b  Z1 


(18) 


(14) 


Therefore,  B  =  0  maps  the  unit  disk  in  the 
with  center  at  Z2  =  -r  and  radius  |  ^  |. 
in  the  Z2  -plane,  if  and  only  if 


z 1  -plane  into  a  disk  in  the  7.^  'Plane 
This  image  will  not  overlap  the  unit  disk 


(20) 


or  equivalently, 

|  a  |  +  |  b  |  <  1 


(21) 
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We  now  try  to  establish  the  same  result  by  using  Theorems  4  and  5.  Making 
the  change  of  variables  according  to  Eqs.  (13),  we  get 

^(Pj*  Vy}  l+P2^"P2"l’Pi  P2  (22) 

H(zr  Z2^  “  F(p1>  P2)  (1+a+b)  +  (1-a+b)  Pj  +  (1+a-b)  p2  +  (l-a-b)  Pj  p^ 

Condition  (2)  of  Theorem  4  dictates  that  F(p^,  1)  =  2(l+a)  +  2  (1-a)  p^  =  0  has 
not  root  in  Re  p^  £  0.  Therefore 

|  a  |  <  1  (23) 

Now  we  tackle  condition  (1)  of  Theorem  4  with  the  help  of  Theorem  5.  We  have 

F(jw,  jft)  =  (1+a+b)  -  (1-a-b)  wft  +  j  £(l-a+b)  w  +  (1+a-b)  S2  J 
and  we  want 


Cq  1  =  (1+a-b)  (1+a+b)  +  (1-a+b)  (1-a-b)  w  >  0, 


(24) 


for  all  real  w 

Let  Aj  =  1+a-b,  A 2  =  1+a+b,  Ag  =  1-a+b,  and  A  4  =  1-a-b.  Then  the  inequality  (24) 
is  valid  if  either  of  the  following  four  conditions  is  satisfied:  (i)  A.  >  0,  A2  >  0, 

Ag  >0,  A^  >  0.  (ii)  A 1  >  0,  A2  >  0,  Ag  <0,  A  4  <  0.  (iii)  A^  <0,  A2  <  0,  Ag  >  0, 

A 4  >  0.  (iv)  A 1  <  0,  Ag  <  0,  Ag  <0,  A4  <  0.  It  is  easy  to  show  that  (i)  is  equivalent 
to  the  inequality  (21),  while  (ii)  -  (iv)  are  incompatible  with  the  inequality  (23). 
b)  A  Special  Class  of  Second -order  Filters 


We  next  consider  the  filter 


H(Zj,  Zg)  l+aZj+bZg+CZ^  Zg 


(25) 
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which  we  shall  call  a  bilinear  filter.  We  first  establish  stability  conditions  in  terms 
of  the  coefficients  a,  b,  and  c,  using  Theorem  3. 

Let 

B(Zj,  Zy)  =  l+azj+bz2+cz  ^  Z2  =  0 

We  have 


1+az  1 

7  2  bfcZj 


(20) 


which  is  a  bilinear  transformation  mapping  circles  into  circles.  The  image  of  the 
unit  circle  8d^  =  (z^‘,  |z^|  =  1)  in  the  z0  -plane  is  then  a  circle.  From  Eq.  (26), 
the  center  of  this  image  circle  is  on  the  real -axis,  and  it  intersects  the  real-axis 
at 

1-a  ,  1+a 

z  =  -  t- -  and  z0  =  -  -t-7- 

2  b-c  2  tH-c 

It  is  easy  to  see  that  condition  (1)  of  Theorem  3  is  satisfied  if 

l&l>‘  <27a> 

and 

1  > 1  <27b> 


Finally,  condition  (2)  of  Theorem  3  is  equivalent  to 
|  a  |  <  1 

The  conditions  (27)  agree  with  those  of  Shanks*. 

We  now  try  to  establish  the  stability  conditions  using  Theorems  4  and  5. 
Making  the  change  of  variables  according  to  Eqs.  (13),  we  get 

If/  „  E(PrP2>  i+p1+p2+p1p2 

(Z1’  Z2*  F(Pl,p2)  B1-H32Pi+B3P2+B4PiP2 


(27c) 


(28) 
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where 


Bj  =  1+a+b+c 
B2  =  1-a+b-c 
B^  =  1+a-b-c 

B4  =  1-a-bfc  (29) 


For  condition  (2)  of  Theorem  4,  we  form 
F(pr  1)  =  2(l+a)  +  2(l-a)Pl  =0 
whence  Re  >  0  means 

|  a  |  <  i 

Now  we  tackle  condition  (1)  of  Theorem  4.  We  have 
F(jw,  jO)  =  Bj  -  B4  wfl  +  j|(B2w  +  B3SI)] 

According  to  Theorem  5,  we  want 

C0, 1  ~  B3  B1  +  B2  B4  W  >0 
for  all  real  w 

which  is  valid  if  either  of  the  following  four  conditions  is  satisfied:  (i)  B^  >  0, 

B3  >0,  B2  >  0,  B4  >  0.  (ii)  Bx  >  0,  B3  >0,  B2<  0,  B4  <0.  (iii)  Bj  <  0,  B3  <0, 
B2  >0,  B4  >  0.  (iv)  Bj  <  0,  B3  <0,  B2  <  0,  B4  <  0. 

It  is  easy  to  show  that  (i)  is  equivalent  to 
|  a+b  |  - 1  <  c  <  1  -  |  a. -b  | 

and  (ii)  -  (iv)  are  ruled  out  by  the  inequality  (30).  It  can  readily  be  verified  that 
(32)  is  equivalent  to  (27). 


(30) 


(31) 


(32) 
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c)  A  Numerical  Example 

Unfortunately,  filters  described  by  Eqs.  (18)  and  (25)  are  the  only  types  for 


which  we  were  able  to  derive  stability  conditions  in  terms  of  the  coefficients  in  the 
z-transforms.  For  more  complicated  filters,  the  application  of  Theorem  3  becomes 
unwieldy.  However,  for  any  filter  with  numerical  coefficients,  we  can  always  apply 
Theorems  4  and  5  and  determine  whether  it  is  stable  or  not  in  a  finite  number  of 
steps.  We  illustrate  by  an  example. 

Consider  the  filter 

TJ/  1  (33) 

H(zl’  z2  ~1  71  71  71  2771  T~ 

1  +  2  Z1  +  2  z2  +  4  Z1  z2  +  4  zl  +  4  Z2 
Making  the  change  of  variables  according  to  Eqs.  (13),  and  proceeding  as 
before,  we  find  that 

FCPj-  !)  =  3p2  +  bpj  +  7 

whose  zeros  have  negative  real  parts.  Therefore,  condition  (2)  of  Theorem  4  is 
satisfied. 

For  condition  (1)  of  Theorem  4,  we  form 
c0,l(w>  4  w  '  1  w  -T 

which  can  certainly  be  negative  for  some  real  values  of  w.  Therefore,  condition  (1) 
is  violated.  The  filter  of  Eq.  (33)  is  unstable. 
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VI.  Concluding  Remarks 

We  have  discussed  in  this  paper  how  to  test  the  stability  of  two-dimensional 
recursive  filters  in  the  frequency  domain.  In  conclusion,  we  would  like  to  point 
out  that  what  we  have  discussed  is  only  part  of  the  stability  problem.  Furthermore, 
it  is  not  the  most  important  part.  The  more  important  questions  we  want  to  ask  are: 
(1)  Hew  do  we  design  filters  that  are  gauranteed  to  be  stable?  and  (2)  If  a  given 
filter  is  unstable,  how  do  we  stabilize  it  without  changing  its  frequency  response? 

One  approach  to  attacking  question  (1)  is  to  study  what  types  of  frequency 
responses  we  can  get  from  classes  of  simple  filters  whose  stability  we  know  how 
to  control. 

2 

With  respect  to  question  (2),  we  might  mention  a  conjecture  of  Shanks  . 

Let  H(Zj,  z2)  =  l/B(z^,  z2)  be  an  unstable  filter.  We  first  determine  a  least  mean- 
square  inverse  of  B(z^,  z-,),  which  we  call 
a  b 

C(z.,  z9)  =  S  _  c  z.m  z9  (34) 

'1  2'  m=0  n=0  mn  12  '  ’ 

where  a  and  b  are  positive  integers. 

Let 

B^,  z2)  z2)  =  1  +  d1Q  z1  +  d  z2  +  dn  zi  z2  . . . 

Then  among  all  polynomials  of  order  (a,  b),  C  is  the  least  mean-square  of  B  if  the 

2  2  2 

coefficients  of  C  are  such  that  (d^  +  d^^  +  d^  +  . . . )  is  minimized.  We  next 

determine  a  least  mean-square  inverse  of  C(z^,  z 2)  which  we  call^z.,  z2).  Then 

Shanks  conjecture  states  that  the  filter  H(Zj,  z 2)  =  l/B^z^,  z 2)  is  stable  and  that 

/\ 

the  magnitude  of  the  frequency  response  of  H  is  approximately  equal  to  that  of  H. 
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The  approximation  becomes  better  when  we  use  larger  values  for  a  and  b  in  Eq.  (34). 

8  9 

It  was  proven  by  Robinson  *  that  a  one-dimensional  version  of  Shanks1  procedure 

does  yield  stable  (one -dimensional)  filters.  However,  whether  the  procedure  yields 

stable  filters  in  two  dimensions  is  still  an  open  question.  Also,  no  analysis  is 

/\ 

available  on  how  close  the  magnitude  of  the  frequency  response  of  H  is  to  that  of  H. 

We  have  yet  another  open  question  with  respect  to  question  (2).  In  the  one¬ 
dimensional  case,  if  we  are  doing  non-real -time  filtering,  then  we  can  always 
decompose  an  unstable  filter  with  no  poles  on  the  unit  circle  into  two  stable  ones 
recursing  in  opposite  directions.  We  associate  the  poles  of  the  original  filter  outside 
the  unit  circle  with  a  filter  recursing  in  the  positive  direction,  and  the  poles  inside 
the  unit  circle  with  a  filter  recursing  in  the  negative  direction.  Is  there  an  analogous 
procedure  in  two  dimension?  It  beats  me. 
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